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ABSTRACT 


The Goncepe of samt cycling Control An power supe irc. 
aad speed ConkLrel Nasmadvanced the apelicavicon Of spheaemuuead 
areas. ON-OFF switching of a thyristor provides a simple and 
economical, Method Ofmeentereo! and resulabaon, This sswitenine 
Peano tLawilsvOor Caulsesethessysten CoO, timig, cycle. 

Basic analysis and design of speed control was performed. 
A describing function was developed to model the power-supply 
and rectifier bridge. Then it was used to predict the fre- 
quency and amplitude of the limit-cycle. 

A digital computer program was used to simulate the sys- 
meomrresponse and to construct the describing function curves. 
femverify the describing function validity, the limit cycle 
Predictions were compared with the simulated results. 

Fourier analyses were performed to determine the ripple 


instability and subharmonic effect of the system output. 
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to LER ODUC RE 


abuabic toni genlavgiens ecugneie | systems have enjoyed considerable 
popularity both jm practice and in’ the Ii cterauureee  nemer ac 
tical reasoms for which limit cycling systems have been 
developed are numerous, not the least important of which is 
wae Simplicity of "am on-off power element. The firsv step 
toward such a theory was actually taken by Lozier (Ref. 6) at 
which time certain qualitative remarks were made. Thereafter 
meereceived Jittle atvention in the@@mglish literature. sie 
paper by Li and Vander Velde (Ref. 7) represented an inde- 
pendent rediscovery of some fundamental principles regarding 
minat cycling control systems., Work of a similar nacure has 
mre, CeCenthy beens published nm the U.S.5.R. 

In Chis arsceusston stress 1s placed on considerations sos 
the limit-cycling system as a control system and its signal 
transmission properties are derived as well as the loop limit 
eve le. 

The discussion is limited to negative Foeduac conten 
rations with a nonlinearity in the ronwerd LOGDRa SUCH aac 
shown in Figure 1.1. Here N stands for the nonlinearity, and 
G(s) for the loop linear elements. It is assumed that this 
is a limit-cycling system. That is, in the absence of an in- 
put, the system sustains a nearly sinusoidal undamped 
@ocillation atv its outpucv. If, as is commonly the case, the 
function of this system is to reproduce the input signal at 


moteone  pONermlevelmabettlS OCluUpUL, UCnene ths TMeosoensc Lome 











PEGURE 2.12 “Sample fimit cy cline centre imeycveme 
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PHCURE 1.2: Typical input and its corresponding output. 
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slow input signal mighs look as Shown Tne hr cure ie ee eis 


well-designed control system the our our hollows Senet a 


the average, within some dynamic Lotlowins errors Overean, 


oenenteyel tank aelavs dhauslat 1, 


Cy. CLers a an approximate model of this 


ore 


error Signal is shown in Ficscure 1535 in tnisemeadeaiaeraee sei 


msoi1d of ampligude 8 a5 vas gocrared with the 1imiteeyele eae 


oe W—-C bias of amplitude @ with the input forced follonim. 


eier Ol’ . 


In this study a system which has properties mentioned 


above will be examined. The specific system to be examined 


is an armature voltage controlled D-C motor. Its armature is 


supplied with a poly-phase, thyristor, A-C to D-C rectifier 


and its field excited by a separate D-C source. 


Twe eoal-an this study as to investigate the control or 


speed of this D-C motor, operating under forced limit eycling 


conditions. A mathematical model is proposed and then ana- 


meacd Using frequency response techniques. 


The thyristor, which is used in the rectifier bridge is 


Boul layer, three 
ivedc is equivalent 


The characteristics 


shown in Figure 1.4. 


two stable modes of 


imtar amount to that 


terminal solid-state device. Functional-— 
ho ltsecas Cube counterpart, a Utmyratwen: 
2a) Cig evinkin sil vere onk Wlavae elaine sticione suas 
in wine torVvardedrrecy mone tie .deyilCemias 
operation, one of high impedance, of sim- 


in the reverse direction and the other, 


of very low impedance. In general the forward break-over 


voltage is inversely proportional to the gate current. This 


iidmeaves than the device could be controlled byecontro Maine 


io 
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FIGURE 1.4: Whyristor transfer characteristics. 
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the level of gate current. The control of the forward break 
over voltage provides the interval, of the forward conducticn 
and thus the amount of power per cycle may be regulated. In 
this paper it is assumed that the entire rectifier bridge is 
eontrolled as a wait rather than by andividual devicesroer 
each thyristor. Wnmevgates of all thyrisvors in thiesrony— 
maese bridge are Connected in parallel and a singe irine 
wena 1S DrOvLdeUmpomo ls Gewrces oom fled “<i remit 
mmaecran 15 contarned in Figure 175. 

Conceptually the system-operaces as Follows 2) sassunic 
that the output speed is below the reference speed. The rec- 
tifier bridge is then supplied with an ON signal, and the 
Brraece output to the armature of the D-C motor is a three= 
paase, full-wave reetified voltage. The filter action of the 
mecvor alternates the harmonics of this voltage and provides a 
melatively smooth speed at the output. When the speed rises 
above the reference level the bridge is supplied with an UF: 
Signal. There is a time delay before the bridge turns off, 
because the shut off conditions for all aye Shores nls (alle 
met. The time delay is approximately of the order of *% 
eycle of the supply frequency. In the OFF state the input 
hor vage to the armature of the D-C motor is zero with reverse 
Bterent essenvially bilockedee te Dbridee wi liieemaim an the 
waves caue Uavil the ©ULpUL Speeduiails below tone reterence, 
whereupon the cycle will repeat itself. This limit cycling 


of the system is the particular area to be considered here. 


ee 
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TI. BASIC ANALYSIS AND DESIGN OF SPEED CONTROL 


D-C motors with ShtW@nt eGxertation are particule sane 
cable for speed control systems.” The speeds ol vite eie vecee mn 


be expressed by 


a Va oe iat 


where, 


Umea Specduol motor 

Va = Armature voltage 

i = Arlagure current 

R = Resistance of armature circuit 
@ = Operating magnetic flux 


Ci —CociiicitenG oF proporvuona lnc. 


It is apparent from equation 2-1, that the speed of such 
movors can be controlled either by varying the armature volit- 
faameia, Of tne operating Tix, © > that is, by Varyane eve 
mee d—-exCitation current of the motor. However, the advan- 
tage of speed control characteristics of a D-C motor cannot 
be fully utilized unless an adjustable D-C voltage can be ap- 
plied to the armature. In order to take full advantage of a 
bee mMOtOr characteristics, it is necessary to provide an ad- 
justable output voltage converter which would convert the A-C 
line voltage into an adjustable D-C voltage, which can be ap- 


piied to the motor. 
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The very popular rotating converter consisting of an A-C 
motor driving a D-C generator, and known as the "WARD LEONARD" 
system, is one of the most widely wsed adjustablle voltage 
Sop veces 

The rectifier converter, combined with complete elec- 
meonwe GCONUCYrOl and, regulating System 1s) pernaps vaecsioat 
practical unit available. Its prime advantage over the ro- 
meine CYpe COnVerters as that only one rovating Machine we 
used instead of three. In addition, the compactness of the 
Emectronic System, the flexibility oe adapvall one yoevarT ens 
SPomalctions, accuracy, and Sensitivity are improved. 

The particular interest here is to control the speed of 
weo-c motor by controlling the conduction angle Of the Ghyrise 
iors, Which are used as an A-C to D-C converter To supply the 
armature of the motor as shown in Figure 1.5. In this config- 
uration a feedback control system was employed to achieve the 
desired performance. 

The system in its simplest block-diagram form is illus- 
feeteda in Nigure 2.1. In order to permit analysis, the syscvem 
components are lumped and partitioned as illustrated in 
ieee 2c, Inethis 1¥elre, block "GS srepresents the morer 
and load. It is assumed that they are entirely linear. 

Block "G,”" is included to take into account any gain that may 
me associated with the comparator, firing circuit, and com-— 
pensation network which may be added to the system. Again, 
this block is assumed to be linear. The source and rectifier 
bridge are represented as block "N" which is assumed to have 


a periodic input and output for development of basic describing 


Ine. 








FiRinG THYRISTOR 
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REFERENCE 
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Pec. ee VOCKNG Patiran Ven ereeilated Speed conerole 





PC a2 one tse Clo OC eC eae ined 
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function (which will be mentioned later in this study). The 
general nature of this block is nonlinear. Block "“H" repre- 
sents a conversion factor of the speed to the voltage. If 
can be considered as a feedback gain. 

The principle of operation can be described as follows: 
mie Speed Sevling yvollage wy, Vas represents the desired motor 
pocece., and the tachemeter fives) Derr data pele ecomm 
Mech 1S directly proportional to the tacho generator specu. 
Pawlebhay On Cf CAcwOMeoueruyOltagGe Is sconiva lent tO pneu rea era 
imeetl,O! Movor Speed.” ine difference between Vp and the 
frehometer yoltage, Krwn, 15 an error signal which can be srep— 
ig@eacnted as 

0 = Vp - KeYm (2-2) 
This signal shows the deviation of the motor speed from its 
desired value. The purpose of the automatic speed control is 
Mommanimize this error Signal and obtain the desired value of 
whe speed. This error Siena l operates the firing Camel: ae 
controls the motor speed by adjusting the conduction angle of 
me thyristors in the rectifier bradge, thus altering the 
Brerare CULrene, andevoltage svalve Of Che arilauurescmncudee 

One procedure for analysis Of The Sysvem 15 by Usineieie 
eenventional transfer function concept. The system transfer 
iinet lon 1s 
Cs) Gi(s)G2(s)N 

Vr(s) cis) Cats ats (2-3) 
The linear portion of the system may be block diagrammed 


as Shown in Figure 2.3. 


af 
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mEGURE 2.3: Block diagram for linear system. 


Tp = ConstTANT 





FIGURE 2.4: _.Motor and load. 
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The transfer function for the mover wanes loacewilte tare 
in Figure 2.6 is developed as follows. The loop equation for 


Pie Soria ville Caliee ier dese 


Va = 1aR ee Yast Kym 


(2-4) 
Torgue-equilibrium equations are 
Tq = Kt ia (2-5) 
From equations 2-5 and 2-6 
me el Nie at co: is 
z (2-7) 
and from equation 2-4 
: eae K alt 
ee eee a ee ee 
L L L (2-8) 


are obtained. 


Choosing state variables as: 


a 


XX} = w 
Xo = 15 
; i oe 
= ge oe Pe oes 
Xa ; X) j X2 
Xo = - Smxy, +Ry, + lv, 
L L L 


A signal flow graph of the system is given in Figure 2.5, 


meeuls 17s indieaving an invegration and X79 and Kan as initial 


fond! F10nNS < 


_ 





EeeGUIR E255: 


Kao Xo 





= Kal 


Signal flow graph for linear system. 
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Using Mason's Gain rude the Cransier function in lap iee. 


MG aio 1 Olea 


Q Mt 
Ge(s) =§ > = IL 
a siete RJ+Lf q+ Rt thE Ky 
JL AUG, (2-9) 
or in Bode form 
te 
Come Rf+tKyKy 
ay 2 RJ+Lf 
——— — 5s Laan SG OF YL 
RE+K,K RE+K, KL “(=ahon 


Numerical values for the motor parameters and load were 


petected from Reference 1, Chapter /, Problem 9. 


They are: 

Armature resistance R = 0.08 ohm 

Armature inductance L = 0.008 Henry 

Motor torgue constant R= 0.95 Newton-meter/Amp 
Motor e.m.f. Sons eane Kae 0.95 Volt-sec/Rad 
Equivalent inertia of movor-load J = 1.38 Ke-m- 

Equivalent VaseOus Uric mon ot f = 0.086 newton-meter 
motor-load ad Sale 
Feedback gain eae 0.318 Volt-sec/rad. 


By using these values and computer program "LISA 360 A PRO- 
GRAM FOR LINEAR SYSTEM ANALYSIS 360 D-16 . 4 . 009", the 
values of magnitude in "db" and phase of the transfer func- 
fron were calculated. From these data Bode plot and Nicmots 
meokewere ObLaimed for Limit cycle investigation. lotus were 


shown in Figures 2.6 and 2.7. 


an 
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Nichols chart for linear part. 
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For a desired value of speed 4l a speci tied 1eadee noni 
practice is to compute the value 70l Faia eae hem eiaieme ei 
produce the desired cutpul, Chenmdestenwche Cirinie ere umnane 
produce a firing angle, %, that biases the rectifier to the 
desired output. Thus, the basic design provides the Gesired 
fesuit in open Loop operation.) Ihe teedback  lcoecpensa sumed 
designed to adjust the system when operating, conditions change. 
If the load’'is changed this causes a change in speed which 
manses an error, and the error is ampiafied and changes i tme 
ifmeeine angie by an amount 8 in apie To counteract the ‘change 
mmespeed. fn such a design it is necessary that the firing 
eeeree lit Sadan, Kae ber Chose sremeenyC bine OCS tbed. SteCagy ote 
speed regulation, then the loop must be checked for stability 
mae: Oor transient response, Since the gain is set on the 
basis of steady state accuracy. It is possible to use an 
open loop approach which includes the nonlinear characteris-— 
ere Of the rectifier. 


a 


Dicom f aoicemwe sen 


Q(s) Gi(s) N Go2(s) 


and 


Q(s) = Ga(s) Vi(s) (22aale) 


Assuming armature voltage is a constant, that is the average 


meswue of the rectifier output, V Using equation 2-3 and 


avé 
motor parameters equation 2-11 becomes 


86 ae 


OS) 2 == ee ROMO Se COST 


a4 








By applying the final value theorem which is 


lim 0(t) = dimers) 
Foo SO 


The steady state value of speed was obtained as below: 


Ww) = 1.045V, 


oo V 


It is known that the average value of the rectifier output is 
emliurnctc10On Of the error sianal, a5 andsune sensitivity or ee 
mmeine Circuit to the error... If the Sensitivity of thesis 
ing circuit is improved by using eatin ay and compensation 
network, better results can be obtained. 

As mentioned before, operation of the system depends on 
existence of a difference between reference speed and measured 
speed. If the "control gain" of the system is defined as the 
tachometer output per volt applied to the amplifier when the 
feedback path is broken. 


Control gain 


ile G om 
Then, 
K, w 
oe cM 
ee ~ Ae 





Thus, the greater the control gain, the closer will be K Win 


CO Ve DOr teh accuracy, Mica Saints required. 
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Til.  PREDICTEON OF LIME ex Chr Cheer l hab ane 


Procedures fOr analysis Ot Cnaemaiticae sporti oneone uae 
system have been previously established by using the trans- 
fer function concept. It is desirable to extend these 
mecnnlques Co the analysis of the nonlinear part oi thes, s— 
mote sAcCGCUAallLY, Sher Tne sieonlinearsperrionetiere Pome ene ra, 
mo exact method for relating the output to the input, in 
erosed form. Most techniques of nonlinear analysis use some 
meer Of approximation, that 1s, assumed linearity of Che 
system over the desired range of operation to be considered. 

The extension of the transfer function representation 
fmmore GifiieulLt.. IG requires thal the nonlinearivy be 
Gescribed by an equation in the frequency domain, and it also 
requires that, this equation be compatible for use with the 
Meamsler functions of the linear components. Practical solu— 
tion to these requirements is the deseribing function technigue 
fee makes Une approximavion in the frequency domaan,. 

The basic philosophy of the describing function is that 
mae Open—loop control system is excited at its inpwt by a 
ame otdal forcing MUunRcEiIOnN. Fhe output cf the nonmlineoarry, 
io periodic and can be represented by its Fourier series. 
fees 1S allustrated in Pigure 3.1. 

Certain assumptions must be made in the process of this 
Megearization, thal the input to the nonlinear component is 
a pure sine wave, and linear portion of the system following 


Mie HOM heAariby 1S a LOW=pass iti terey which 1s Common ty aeruc 
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in most physical feedback control systems thse ees eoccumes 
that the only terms of Che Sowrier soem eo ia vee Gomes mrt. 
the output of the system are the D-C term and the fundamental 
hrequency term. The fundamental frequency, term 1 thew sourrer 
series has the same frequency as the input signal but may dif- 
mer in amplitude andeohase. When the “oop is closedrumsc 
harmonic will propagate itself around the loop and result in 
ae continuous oscillation of the system. Then the describing 
function may be defined as the ratio of the fundamental term 
in the Fourier series for the output wave to the input sinu- 
gonad considering bouh magnitude and phase. 

lne-order to @eu the describing samc ti On Lor auynonla head 
ie che OUbpDUL of this nonlinearity as represented by ies 


Fourier series. 


A 
f(t) = s + A,Coswt + BySinwt + --- + A Cosnut + B_Sinnot 
(3-1) 
The coefficients of the series are given by 
5 a 
a Fe 1G) OO eia0 ic ili (3-2a) 
O 
5 ah 
BL ar fe yee oa WilaKine Cle (3-2b) 
O 


temecciahavion OF Eme escribane Tuncolon requires Thav A, 
and B for n ? 2 are negligible. Thus the coefficients of 
interest are A, and Bi, and they define the fundamental fre- 


avency OUTDUL as 


O(t) = | A,? +4 B,2|? /_tan™ (A1/B1) C= 


eA (| 








Then the describing function for tiers ys vee aoe tt we ie Lew 
De Pa eure. 3 ales 
2 Za 
[Ay + Ba | 


oS eee /tan '(A1/B1) - ¢ . 
B (3-6) 


The above criteria will now be applied to the system 
moder consideration. 

A mathematical description of the output wave shape of 
mae nonlinear block must first be obtained for dérivation of 
the describing function. After that the Fourier analysis is 
forked out to get sinusoidal components of it. 

hReterring to Pigure 3.2, the ouwlout of the system hacwa 
eerbcain constant speed, and @ ripple which is superimposed on 
mars Constant level. The first harmonic of the ripple is de— 
fined as BSinwt, and the difference between the constant level 


and the reference speed as 


The system is presumed to operate as follows: 

a) When the constant level of the speed plus the ripple 
is less than the reference level the rectifier bridge is "ON" 
between time a to b, and the output of the nonlinear portion 
of the system is a three-phase, full-wave rectified voltage 
@imimean Value A. 

b) At time b where the constant speed plus the ripple 
equals the reference speed, the rectifier bridge receives an 


BOorH signal . 


28 





BSinwt 







f(t) =zAnCosnwt+ZBnSinnwt 


PIG ES. 


1 : 
I ! a 
( 
| ee iam 
| 
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ec) The rectifier is turned off at time equal to c, a 
delay of 1/240 second after time b. 

d) The rectifier remains off until the output speed 
again falls below the reference. (Time c to T.) ‘Then the 


output of the nonlinear block may be summarized as follows: 





Pr 0 O0<t<a 
= A ass) 
- c-t 
ey aE D=t<c 
= 0 e<ccl 


The times a, b, and c can be defined in terms of the 


minst harmonic of the ripple. 


ears -* 
= 
Ww 
TT 
a=--—-— -)D 
= Gd 
hha 


To find the coefficients of the Fourier series equations 
Secea,0 are used. The integration must be evaluated over one 
full cycle of the disturbance frequency, T. 


Prece@dure is as Tollows: 


B ae 
A. = P(t) Cocnur dt 
O 
c-t _ 
—— = 240b + 1- 2uOo t 
e-b 
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Cc c CG 
[ Coswt ate2uon [ Coswt ac-ato [ tC Ces on dt 
a b 


imcesratings and .subs taGurine Suric sums: 


40 


2 che ra omcnine earns a eaeiie 


es a 





—-24O0cSinwet 2 Coswb+240bSinwb 





Cancelling terms, the above reduces to: 


240 


A. = eA - Sinwbt = (Coswb-Coswe ) 


1” Tw 








mews trigonomatric 1dentities and substituting the follow 


ing values, 


Bie eo 
= = = and Coswb = 
ee 5, Sinwb = RB an Oswb = B 
The final form of the A, term is found as equation 3-5. 
a A oy Vv 2 2» G) r Od 
A, = 5 —-O + Fi | B“-a 3-008 2 Ooi D0 | 
(3-5) 
Proceeding in a like manner B, is found as equation 3-6. 
eee ees pee) 2 40 W : W 
By =) B = a a c [1-cos dO 0 | vi - a2 ai DO 0 
(3-6) 


Bel 





The describing Tunct longo elem. 


en 2 2 
1D 
iti a oe 


B (3-7) 


where 


Pane, (i. /Ba) 


cD 
i 


Equations 3-5, 3-6, and 3-7 are suitable for evaluating 
mae numerical value of the describing function. Due to their 
complex nature, the digital computer was used to overcome the 
labor required to calculate more than a limited number of 
points. 

The computer program with numerical data is contained in 
mopendix A. 

Since the describing function, which is derived above is 
dependent upon the frequency, the magnitude of the input, B, 
and the difference between constant level of the output and 
ieee lCrence Specd Op. separate sets Of calculatwens are 
meeessary for each choice of an independent variable. Hor 
eemvenience the magnitude of the ianput was chosen as the ins 
Mependent variable. The flow chart for the computer program 
meeshown in Figure 3.3. 

Tommake tne —Ppredice lO women tweyvele —Characverisuac  vmc 
describing functions were drawn on the Nichols chart. They 
are shown in Figures 3.4 through 3.6. 

The values of alpha and frequencies for evaluation of 
describing functions were obtained from the simulation of the 


system which will be discussed later. 
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B= 8B+DeELTAB 








HAVE REQUIRED 
NUMBER of PoinTsS BEEN 
Computed ? 


PRINT OvuTPUT 


B,WN,-!/N, 








ANGLE. 


FIGURE 3.3: Flow chart to compute value of describing 
Fue tron. 
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mIGURE 3.4: Describing Eunetaon (-1/N} for Limit cycling 
thyristor rectifier bridge. 
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PEGURE 3.5: Describing stunctiom {(~-l/N) for limit cycling 
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These values were tabulated in Table 3-l. 


Vv REFRRENCE MEAN VALUE 
i Om IN RAD/SEC OF OUTPUT ALPHA FREQUENCY 


31 930 rpm 97.4882 97.4843 0.0055 56 .61HZ 
32.8% 984 rpm ores 103.1378 0.0065 50.15HZ 
33.5. 1005 rpm 105.3459 105.337 0.0089 Override 


TABLE 3-1. Tabulation of simulation results for ve equal 
Sloe omens. > a ollus. 


Having derived the vransfer function for the linear per 
tron Of the system and describing function for the nonlinear 
part of the system there are now a number of ways to deter- 
ime the condition necessary for limit cycling. 

The characteristic equation for the system is written 


from equation 2-3 as: 


alee ote Gi(s) Gos) iss) N = QO 


This can be rearranged as: 


G,(s)Go(s)H(s) = - ~ = G(s) 


Thus the intersection of a graphical representation of G with 
-1/N determines the critical point at which limit cycling will 
meeur. Either a Nyquist or a Nichols plot is convenient for 
Um purpose. Ine choice ot the Nicholls pilovt tor these ana 


AySes 15 arditrary. 


oT 








To determine the docation lo! time init ee yelearnemde. 
ecribing function curves were superimposed on the frequency 
response curve of the linear part for frequency 39.1/7 Hz 
with alpha 0.0089, for frequency 50.15 HZ with alpha 0.0065, 
and for frequency 56.61 HZ with alpha 0.0055. The graphs 
ere illustrated in Figures 3.7 through 3.9. 

The values of frequency and amplitude were measured at 


intersections of G@ with -1/N. The results are tabulated in 


Table 3-2. 
Vv Alip iigiem el  sticeed lemMevwesot Freamanen of 
r -1/N -1/N G -: 
Gel 0.0063 Oona Enz 56.65 HZ 
BO O50072 Dom Hz 50m 1s ga 
bo! 0.0095 39.17 HZ 39.23 Ha 


femies 35-2. Tabulation of the ebserved limit cycle frequency 
and amplitude. 

From the comparison of the results in Table 3-2, it can 
be seen that the frequency of the -1/N curve and frequency of 
mae CULrve are very lose CoOveach Ovner this is true for 
pele three cases. 

These results were obtained when alpha was equal to the 
melues thay are indicaved in Table 3-1. Corresponding Go 
these alpha values the amplitude of the limit cycle was found 


as indicated in Table 3-2. 
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Thus there is good agreement with simulation results in 
both frequencies of the limit cycle an@ the values of ealpha 
predicted by the describing function. This indicates the 
Validity of the Gdeseribinge functions a predict ionmoms ule 


Limit eycle of the system. 


Ho 








IV. SIMULATION OF THE SYSTEM AND SUBHARMONIC ERPECTS 


To determine time vada yey ac eclers creat iia cumme creme 
and Go See the perProriianceme. ime wens Colmer me time domain, 
the system was simulated on the digital computer. The flow 
emart for the simulatiom is shown in Ficure 4.1. “Naval 
Postgraduate School computer library subroutine "INTEG 2" 
was used to solve the differential equations required in the 
Simulation. The complete program is contained in Appendix B. 

the supply volvage was Considered to be a three-phase 
60 cycle A-C source with a peak amplitude of 115 volts through- 
eur the Simulation. 

The forcing voltace for the system, that 18 applied vo 
the armature of the D-C motor is a three-phase, full-wave, 
mectitied voltage during time @ to b, a@ portion of a sine 
wave for the appropriate conduction phase during time b to 
ec, and zero voltage with no reverse current through the in- 
mcwance during c to TF of Figure 3.2. 

Graphical output from this program are shown in Figure 4.2 
imebiictrabing the input voltage wave form, and Figures 473 
through 4.7 illustrating the output speed wave forms. 

Simulation was performed by changing the reference from 
930 r.p.m. to 1260 r.p.m. These are shown in Table 4-1. 

When a linear system is driven with a sinusoidal input 
it is known that the steady state output is also a sinusoid, 
differing from the input in amplitude and phase but having 


precisely the same frequency. For linear systems, that is 


nS 








a CALL 
eee 






Form B3-PHASE Fur.-~WAVE 


RECTIFIED VOLTAGE 


DETERMINE ConwpucTiONn 


PHASE A,B,C. 


ERROR = Out -~ REFERENCE 









NO Is £RROR Positive ? YES 


WAS Vin For LAS 
NO YES 
ITERATION EQUAL 


TO ZERO? 





Vin = A Ximun 
EHOSE Hoot. 






Vinz ConDuctTinG V., = ZERO 
Purns=e A,B,C. 


SO_rve Digge. 
EQUATIONS 


SOLVE Dire. 
EQUATIONS 


SOLVE DIFE. 
EQUATIONS 







GO TO 4 
CON DUCTION 


PHASE VOLTAGE EQUALS 
ZERO 


| vas 00 
Vin = ZERO 


we 


FIGURE 4.1. Flow chart for simulation. 
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Simulation results, output speed. 
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FIGURE 4-4. Simulation results, output speed 
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FIGURE 4-5. Simulation results, output speed. 
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FIGURE 4=6. Simulation results, outputs speed. 
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DOES SUB- 


REFERENCE PER POD LOE FREQUENCY OF HARMONIC 
VOLTAGE SPEED iN ieee ls LINED CYCLE EAL Sik 
31 volt 930 rpm 0.017663 Sec 56. 6 lueelZ ES 
32 fewer 960 rpm .017631 Sec SiO. (ae Yes 
Be. 3 volt 969 rpm -O021173 Sec 47, QO Z. YES 
Boe) VOLT 975 rpm SOLo9 24 ose Oe eee YES 
B220 Volt 978 rpm nO 5 sec Be o), He vis 
32.8 volt 984 rpm 0.01994 Sec 50.13 HZ YES 
Bee Volt 990 rpm "Oescil” Sec Hi, O02 HZ ‘ins 
Be. 3 VOLt 999 rpm nO 3S Sec Hi 76 We YES 
Bo. volt. 6 ©1005 rom "025525 Sec Sorel (fae YES 
Boeo volt 1014 rpm .02545 Sec 39.293 HZ Yes 
fee volt 1020 rom 1027 036m Sec 5 ore 3 HZ ies 
Boe volte 1050 rom .0394 Sec 25.4 Le Nene 
ey. Volt 1110 rpm NO Peer CYCLE 
wee volt 1200 rpm NO ete eae CYCLE 
Wee VOLG 1260mrom NO LIMIT CYCLE 
HABE l- Tabulationeor caimmliation results. 
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true at all frequencies. In the case of a nonlinear system 

a sinusoidal input does not guarantee an output which is a 
sinusoid. In the usual case, the output wave contains har- 
monic frequency components which are integral multiples of 
the forcing frequency. In some cases, the output wave may 
contain one or more frequencies which are lower than the for- 
cing frequency. Characteristically, these frequencies are 
integral submultiples of the forcing frequency which are 
called subharmonic frequencies. The amplitude of the sub- 
Harmonic component may be small compared with the amplitude 
of the component at the driving frequency, or it may be so 
faree that the driving frequency may be neglected for” partic— 
ular computations. 

If the output wave shapes in Figure 4.3 through 4.7 are 
inspected, it is observed that the wave forms have a ripple 
instability which is smail in amplitude compared to the ave- 
rage ripple amplitude. One can think this Seay is 
ferusea by SUbDharmonic eifects of the system. It can be ob= 
servea by inspection, but it 1S not possible to obtain its 
frequency and its amplitude from the output wave shapes, be- 
cause they do not have a regular pattern. To figure out this 
miacability, the Pourier andlysis of the output wave shapes 
for reference speed 930, 984, and 1005 r.p.m. were made. 
During the Fourier analysis, several limit cycle periods were 
taken as a fundamental period of the Fourier, then the same 
harmonics of these fundamental frequencies which are also sub- 
HarntOnLes Of rhe lime cycle were obvained, | rom Tocse swe 


sults, amplitude and the order of subharmonic components of 


De 








the limit .cycle in the output wave were calculated. These 
were tabulated in Table 4-2 through 4-4. 

In Figure 4.8 through 4.10, the amplitude of the sub- 
harmonics was drawn to see the relative heights of the com- 
ponents. In Figure 4.8, it is seen that, at reference speed 
930 r.p.m. 1/8 and 1/16 subharmonics have larger amplitude 
compared to others. Their amplitudes are approximately 15.6 
percent of the output ripple amplitude. In Figure 4.9 which 
was drawn for reference speed 984 r.p.m., 1/2 subharmonic 
fas larger amplitude and it is equal 25.3 percent of the out— 
put ripple. In Figure 4.10, which is for reference speed 
mOO5 r.p.m. the largest subharmonic is 1/7. It has amplitude 


16.1 percent of the ripple amplitude. 
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V. CONCLUSION ANDSRE COMMENDATION DRORs aCe ee ba 


The main conclusion that might be reached in reviewing 
the experimental results is that the D-C motor at constant 
iecld current operates au an almosy scons uaneacp ec dam idcmmanre 
speed may be controlled over wide ranges by armature-voltage 
@omtrol, This can be seen from Pigures 4.3 through 4.7. “i 
moese figures are inspected closely it can be observed that 
the ripple amplitudes are very small compared to the output 
Value. They are approximately within the + 0.025 percent of 
the desired value of the output. This may be acceptable for 
Darel vie Appl calaona an une iInduscry - 

Armature-voltage control was performed by adjusting the 
MmoNdues Tom phase ol the thyristvor-recviiier bridge. waais 
eeuises the limit cycling in the control loop. 

Tne descrubing funerion cerivea in Chapter Tift iereevEdes 
a mean by which the nonlinearity of the system can be repre- 
sented and the overall performance of the system analyzed. 
Mees POSsSible to predict the approximate frequency and ampii-— 
Mrcic Of the limit cycle as illustrated in Figures 3./ through 
ieee in Order to predict the response more precisely, 10 as 
necessary to predetermine the value of alpha, that might be 
Aone by simulating the system and measuring the output mean 
value then taking the difference between output mean value and 
reference. 

One tentative conclusion that can be achieved by examina- 


wren Ol tie experimental data is that, under certain conditions 


60 








the forced-limit-cycling speed control system will have sub- 
harmonie ripple instabi tii. 

The describing function derived aecounts  semiy stemware 
first harmonie of the Ouppuy vollagem rome pic enikieeer Sas 
tion of the system and assumes that higher harmonics are 
negligible. However, as seen in Figures 4.3 through 4.7 the 
euher harmonics can exist. in some cases these harmonics 
ean be submultiples of the limit cycle frequency. 

In this study Fourier analysis was performed to see 
whether the subharmonic exists or not. As a result Tables 
W-2 through 4-4 and Figures 4.8 through 4.10 were obtained. 
These show that the ripple instability is caused by subhar- 
meme Cifects of the system nonlinearity. But this meme 
good enough to predict the existence of the subharmonic in 
ime OUTDUL Of the system. It can five some insight about 
the nature of the problem. There are no general rules avail- 
mpee Whach Geline the conditions necessary for the occurrence 
of the subharmonic oscillation. For mathematical analysis 
several of the references noted in the hist of Reierences 
suggest that a possible method of prediction of the subhar- 
monic ripple instability in the system would be to employ a 
Tie inoue Describing Function. Any further Sscudy should 
give some consideration to these analytical methods to under- 
evuand ripple instability and existence of the subharmonic. 

Another area for study is to increase the gain and see 
the effect of the gain on the ripple instability. This may 
require some sort of compensation network to eliminate the 


Pisvaoikivey im une OULDUL for a range ol gain values. 
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APPENDIX A 


Computer program for computing value of describing 
curves. Numerical data for -1/N curves plotted in Figure 3.4 
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